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Abstract 

Wireless networks can be self-sustaining by harvesting energy from ambient radio-frequency (RF) signals. 
Recently, researchers have made progress on designing efficient circuits and devices for RF energy harvesting 
suitable for low-power wireless applications. Motivated by this and building upon the classic cognitive radio (CR) 
network model, this paper proposes a novel method for wireless networks coexisting where low-power mobiles in a 
secondary network, called secondary transmitters (STs), harvest ambient RF energy from transmissions by nearby 
active transmitters in a primary network, called primary transmitters (PTs), while opportunistically accessing the 
spectrum licensed to the primary network. We consider a stochastic-geometry model in which PTs and STs are 
distributed as independent homogeneous Poisson point processes (HPPPs) and communicate with their intended 
receivers at fixed distances. Each PT is associated with a guard zone to protect its intended receiver from ST's 
interference, and at the same time delivers RF energy to STs located in its harvesting zone. Based on the proposed 
model, we analyze the transmission probability of STs and the resulting spatial throughput of the secondary network. 
The optimal transmission power and density of STs are derived for maximizing the secondary network throughput 
under the given outage-probability constraints in the two coexisting networks, which reveal key insights to the 
optimal network design. Finally, we show that our analytical result can be generally applied to a non-CR setup, 
where distributed wireless power chargers are deployed to power coexisting wireless transmitters in a sensor network. 

Index Terms 

Cognitive radio, energy harvesting, opportunistic spectrum access, wireless power transfer, stochastic geometry. 

I. Introduction 

Pov^ering mobile devices by harvesting energy from ambient sources such as solar, w^ind, and kinetic activities 

makes wireless networks not only environmentally friendly but also self-sustaining. Particularly, it has been reported 

in the recent literature that harvesting energy from ambient radio-frequency (RF) signals can power a network of 

low-power devices such as wireless sensors [l]-[6]. In theory, the maximum power available for RF energy 

harvesting at a free-space distance of 40 meters is known to be 7uW and luW for 2.4GHz and 900MHz frequency. 
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respectively [1]. Most recently, Zungem et al. have achieved harvested power of 3.5mW at a distance of 0.6 meter 
and luW at a distance of 11 meters using Powercast RF energy-harvester operating at 915MHz [1]. It is expected 
that more advanced technologies for RF energy harvesting will be available in the near future due to e.g. the rapid 
advancement in designing highly efficient rectifying antennas [3]. 

In this work, we investigate the impact of RF energy harvesting on the newly emerging cognitive radio (CR) 
type of networks. To this end, we propose a novel method for wireless networks coexisting where transmitters 
from a secondary network, called secondary transmitters (STs), either opportunistically harvest RF energy from 
transmissions by nearby transmitters from a primary network, or transmit signals if these primary transmitters 
(PTs) are sufficiently far away. STs store harvested energy in rechargeable batteries with finite capacity and apply 
the available energy for subsequent transmissions when batteries are fully charged. The throughput of the secondary 
network is analyzed based on a stochastic-geometry model, where the PTs and STs are distributed according to 
independent homogeneous Poisson point processes (HPPPs). In this model, each PT is assumed to randomly access 
the spectrum with a given probability and each active (transmitting) PT is centered at a guard zone as well as a 
harvesting zone that is inside the guard zone. As a result, each ST harvests energy if it lies in the harvesting zone 
of any active PT, or transmits if it is outside the guard zones of all active PTs, or is idle otherwise. This model is 
applied to maximize the spatial throughput of the secondary network by optimizing key parameters including the 
ST transmit power and density subject to given PT transmit power and density, guard/harvesting zone radius, and 
outage-probability constraints in both the primary and secondary networks. 

Our work is motivated by a joint investigation of the proposed conventional opportunistic spectrum access and 
the newly introduced opportunistic energy harvesting in CR networks, i.e., during the idle time of STs due to the 
presence of nearby active PTs, they can take such an opportunity to harvest significant RF energy from primary 
transmissions. Specifically, as shown in Fig. 1, each ST can be in one of the following three modes at any given 
time: harvesting mode if it is inside the harvesting zone of an active PT and not fully charged; transmitting mode 
if it is fully charged and outside the guard zone of all active PTs; and idle mode if it is fully charged but inside 
any of the guard zones, or neither fully charged nor inside any of the harvesting zones. 

A. Related Work 

Recently, wireless communication powered by energy harvesting has emerged to be a new and active research 
area. However, due to energy harvesting, existing transmission algorithms for conventional wireless systems with 
constant power supplies (e.g., batteries) need to be redesigned to account for the new challenges such as random 
energy arrivals. For point-to-point wireless systems powered by energy harvesting, the optimal power-allocation 
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Fig. 1. A wireless energy harvesting CR network in which PTs and STs are distributed as independent HPPPs. Each PT/ST has its intended 
information receiver at fixed distances (not shown in the figure for brevity). ST harvests energy from a nearby PT if it is inside its harvesting 
zone. To protect the primary transmissions, ST inside a guard zone is prohibited from transmission. 

algorithms have been designed and shown to follow modified water- filling by Ho and Zhang [7] and Ozel et al. [8]. 
From a network perspective, Huang investigated the throughput of a mobile ad-hoc network (MANET) powered 
by energy harvesting where the network spatial throughput is maximized by optimizing the transmit power level 
under an outage constraint [9]. 

Among other energy scavenging sources such as solar and wind, background RF signals can be a viable new 
source for wireless energy harvesting [10]. A new research trend on wireless power transfer is to integrate this 
technology with wireless communication. In [1 1] and [12], simultaneous wireless power and information transfer has 
been investigated, aiming at maximizing information rate and transferred power over single-antenna additive white 
Gaussian noise (AWGN) channels. For broadcast channels, Zhang and Ho have studied multi-antenna transmission 
for simultaneous wireless information and power transfer with practical receiver designs such as time switching and 
power splitting [13]. Moreover, Zhou et al have proposed a new receiver design for enabling wireless information 
and power transmission at the same time, by judiciously integrating conventional information and energy receivers 
[14]. For point-to-point wireless systems, Liu et al. have studied "opportunistic" RF energy harvesting where the 
receiver opportunistically harvests RF energy or decodes information subject to time-varying co-channel interference 
[15]. More recently, Huang and Lau have proposed a new cellular network architecture consisting of power beacons 
deployed to deliver wireless energy to mobile terminals and characterized the trade-off between the power-beacon 
density and cellular network spatial throughput [16]. 
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In another track, the emerging CR technology enables efficient spectrum usage by allowing a secondary network 
to share the spectrum licensed to a primary network without significantly degrading its performance [17]. Besides 
active development of algorithms for opportunistic transmissions by secondary users (see e.g. [18], [19] and 
references therein), notable research has been pursued on characterizing the throughput of coexisting wireless 
networks based on the tool of stochastic geometry. For example, the capacity trade-offs between two or more 
coexisting networks sharing a common spectrum have been studied in [20]-[22]. Moreover, the outage probability 
of a Poisson-distributed CR network with guard zones has been analyzed by Lee and Haenggi [23], where the 
secondary users opportunistically access the primary users' channel only when they are not inside any of the guard 
zones. 

B. Summary and Organization 

In this paper, we consider a CR network with time slotted transmissions and PT/ST locations modeled by 
independent HPPPs. The ST transmission power is assumed to be sufficiently small to meet the low-power 
requirement with RF energy harvesting. Under this setup, the main results of this paper are summarized as follows: 

1) We propose a new CR network architecture where STs are powered by harvesting RF energy from active 
primary transmissions. We study the ST transmission probability as a function of ST transmit power in the 
presence of both guard zones and harvesting zones based on a Markov chain model. For the cases of single- 
slot and double-slot charging, we obtain the expressions of the exact ST transmission probability, while for 
the general case of multi-slot charging with more than two slots, we obtain the upper and lower bounds on 
the ST transmission probability. 

2) With the result of ST transmission probability, we derive the outage probabilities of coexisting primary and 
secondary networks subject to their mutual interferences, based on stochastic geometry and a simplified 
assumption on the HPPP of transmitting STs with an effective density equal to the product of the ST trans- 
mission probability and the ST density. Furthermore, we maximize the spatial throughput of the secondary 
network under given outage constraints for the coexisting networks by jointly optimizing the ST transmission 
power and density, and obtain simple closed-form expressions of the optimal solution. 

3) Furthermore, we show that our analytical result can be generally applied to even non-CR setups, where dis- 
tributed wireless power chargers (WPCs) are deployed to power coexisting wireless information transmitters 
(WITs) in a sensor network, as shown in Fig. 2. Practically, WPCs can be implemented as e.g. wireless 
charging vehicles [24], or fixed power beacons [16] randomly deployed in a wireless sensor network. Based 
on our result for the CR network setup, we derive the maximum network throughput of such wireless powered 



5 



O 








® 




o0 










A 



A WPC: Transmitting 
O WIT: Harvesting 
© WIT: Transmitting 

O WIT Idle 

(not fully charged) 
O Harvesting zone 



o 




Fig. 2. A wireless powered sensor network in which WPCs and WITs are distributed as independent HPPPs. Each WIT has intended 
receiver at a fixed distance (not shown in the figure for brevity). WIT harvests energy from a nearby WPC if inside its harvesting zone. 
Unlike the CR setup in Fig. 1, the guard zone is not applicable in this case, and thus a fully charged WIT can transmit at any time. 



sensor networks in terms of the optimal density and transmit power of WITs. 

The remainder of this paper is organized as follows. Section II describes the system model and performance 
metric. Section III analyzes the transmission probability of energy-harvesting STs. Section IV studies the outage 
probabilities in the primary and secondary networks. Section V investigates the maximization of the secondary 
network throughput subject to the primary and secondary outage probability constraints. Section VI extends the 
result to the wireless powered sensor network setup. Finally, Section VII concludes the paper. 

Notations of selected symbols used in this paper are summarized in Table I. 



II. System Model 



A. Network Model 



As shown in Fig. 1, we consider a CR network in which PTs and STs are distributed as independent HPPPs 
with density and A^, respectively, with A^ <C A^. It is assumed that time is slotted and each PT independently 
accesses the spectrum with probability p at each time slot. Thus, the point process of active PTs forms another HPPP 
with density A^ = pXp, according to the Coloring Theorem [27], which varies independently over different slots. 
For convenience, we refer to active PTs simply as PTs in the rest of this paper. We denote the point processes 
of PTs and STs as $p = {X} and $s = {y}, respectively, where X, y G denote the coordinates of the 
PTs and STs, respectively. In addition, it is assumed that each PT/ST transmits with fixed power to its intended 
primary/secondary receiver (PR/SR) at distances dp and dg, respectively, in random directions. We denote the fixed 
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TABLE I 
Symbol Notation 



Symbol 


Meaning 




T TT~)T~)1~) ^-P TIT" CT" 

Hrrr or rl, 1^1 


Ap, As 


Density oi PT, ST 




Point process oi transmitting ST 


ZD ZD 


Transmission power of PT, ST 


a 


Path-loss exponent 


2 


AWGN noise power 


^^.^/^ 


Radius of guard zone, harvesting zone 




Union of guard zones, harvesting zones 


Pg^Ph 


Probability of being outside guard zone, inside harvesting zone 




Maximum number of energy-harvesting time slots to fully charge the battery 


b[I ,x) 


Disk with radius x > centered at T G 




Annulus with radii < x < y centered at T G M"^ 


H 


Distance between a typical ST and its nearest PT 


Pf^Pt 


Probability of a ST being fully-charged, transmitting 


X , y 


lypical rl, 51 




Typical PT, ST corresponding to the respective PR/SR at the origin 


Op, Os 


SINR threshold for the primary network, secondary network 


p(p) p(s) 

-'out 5 -'out 


Outage probability for the primary network, secondary network 


Cs 


Outage probability constraint for the primary network, secondary network 




Spatial throughput of the secondary network 



transmission power levels of PTs and STs as Pp and Pg, respectively. We assume Pp » Ps in this paper for energy 
harvesting applications of practical interest. 

STs access the spectrum of the primary network and thus their transmissions potentially interfere with PRs. To 
protect the primary transmissions, STs are prevented from transmitting when they lie in any of the guard zones, 
modeled as disks with a fixed radius centered at each PT. Specifically, let 6(r, x) C M? represent a disk of radius 
X centered at T G M?; then b{X,rg) denotes the guard zone with radius rg for protecting PT X G Define 
G = Uxg^p K^^'^g) union of all PTs' guard zones; accordingly, a ST y G $s cannot transmit if Y e Q. 

Note that in practice the guard zone is usually centered at a PR rather than a PT as we have assumed, while our 
assumption is made to simplify our analysis, similarly as in [18]. We further assume dp <^ Vg to guarantee that 
guard zones centered at PTs (rather than PRs) will protect the primary transmissions properly. Under the above 
assumptions, the probability pg that a typical ST, denoted by Y^, does not lie in Q is equal to the probability that 
there is no PT inside the disk centered at Y^ with radius r^, i.e., b{Y^^rg). Note that the number of PTs inside 
b{Y~^^ Tg), denoted by N, is a Poisson random variable with mean nrgXp; thus, its probability mass function (PMF) 
is given by 

Py{N = n} = e-^^g^^ ^ , n = 0, 1, 2, ... (1) 

n! 
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Consequently, pg can be obtained as 



(2) 



= PrjAT = 0} 



(3) 



(4) 



We assume flat-fading channels with path-loss and Rayleigh fading; hence, the channel gains are modeled as 
exponential random variables. As a result, in a particular time slot, the signals transmitted from a PT/ST are 
received at the origin with power gxPp\X\~^ and gYPs\Y\~^^ respectively, where {gx}xe^p and {^y}yG<i>. are 
independent and identically distributed (i.i.d.) exponential random variables with unit mean, a > 2 is the path-loss 
exponent, and |y| denote the distances from node X^Y to the origin, respectively. 

B. Energy-Harvesting Model 

To make use of the RF energy as an energy-harvesting source, each RF energy harvester in a ST must be 
equipped with a power conversion circuit that can extract DC power from the received electromagnetic waves 
[2]. Such circuits in practice have certain sensitivity requirements, i.e., the input power needs to be larger than 
a predesigned threshold for the circuit to harvest RF energy efficiently. This fact thus motivates us to define the 
harvesting zone, which is a disk with radius centered at each PT X G $p with ^ Tg. The radius is 
determined by the energy harvesting circuit sensitivity for a given P^, such that only STs inside a harvesting zone 
can receive power larger than the energy harvesting threshold, which is given by Ppr'^^. The power received by 
a ST outside any harvesting zone is too small to activate the energy harvesting circuit, and thus is assumed to be 
negligible in this paper. 

Let Vh) represent the harvesting zone centered at PT X G such that a ST F can harvest energy from 
one or more PTs if Y e H, where H = [Jxe^p ^h) denotes the union of the harvesting zones of all PTs. 
The probability ph that a typical ST Y^ lies in T-L is equal to the probability that there is at least one PT inside the 
disk h{Y~^^rh). Similar to (1), the number of PTs inside h{Y~^^rh), denoted by K, is a Poisson random variable 
with mean vrr^Ap and PMF given by 



Pr{i^ = k} = e 




A: = 0,l,2,... 



(5) 



k\ 
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Accordingly, is given by 



Ph = Pr{y* G n} 



(6) 



= Pr{K > 1} 



(7) 



OO 




= 1 _ g-^rr^^p 



A;! 



(8) 



(9) 



Since Xp and r/^ are both practically small, we can assume vrr^Ap ^ 1. Thus, ph given in (8) can be approximated 
as Pr{K = 1} by ignoring the higher-order terms with k > 1. Therefore, when G is inside the harvesting 

zone of one single PT most probably, which equivalently means that the harvesting zones of different PTs do not 
overlap at most time. As a result, the amount of average power harvested by G ?^ in a time slot can be 
lower-bounded by r]PpR~^ where R < denotes the distance between Y^ and its nearest PT, and < 77 < 1 
denotes the harvesting efficiency. Note that the harvested power has been averaged over the channel short-term 
fading within a slot. 

C. ST Transmission Model 

We assume that each ST has a battery of finite capacity equal to the minimum energy required for one-slot 
transmission with power Ps for simplicity. Upon the battery being fully charged, a ST will transmit with all stored 
energy in the next slot if it is outside all the guard zones. We denote the probability that Y^ has been fully charged 
at the beginning of a time slot as and the probability that it will be able to transmit in this slot as pt. As 
mentioned above, the point process of PTs $p varies independently over different slots, and thus the events that 
a ST has been fully charged in one slot and that it is outside all the guard zones in the next slot are independent. 
Consequently, pt can simply be obtained as 



where pg is given in (4), and pf will be derived in Section III. 
D. Performance Metric 

For both PRs and SRs, the received signal-to-interference-plus-noise ratio (SINR) is required to exceed a given 
target for reliable transmission. Let Op and Og be the target SINR for the PR and SR, respectively. The outage 
probability is then defined as P^^J, = Pr{SINR(^) < 9p} for the primary network and P^^t = Pr{SINR(^) < 6>J 

(v) (s) 

for the secondary network. The outage-probability constraints are applied such that P^ut ^ ^out ^ with 

given < 6p, 65 < 1. Note that the transmitting STs in general do not form a HPPP due to the presence of guard 



Pt = PfPg 



(10) 



zones and energy harvesting zones, but their average density over the network is given by pt^g. Accordingly, given 
fixed PT density Xp and transmission power Pp, the performance metric of the secondary network is the spatial 
throughput Cs (bps/Hz/unit-area) given by 

Cs=PtXslog2{l + 9s), (11) 

under the given primary/secondary outage probability constraints and e^. 

III. Transmission Probability of Secondary Transmitters 

In this section, the transmission probability of a typical ST pt given in (10) is analyzed using the Markov chain 
model. For convenience, we define M as the maximum number of energy-harvesting time slots required to fully 
charge the battery of a ST. Since the minimum power harvested by a ST in one slot is rjPpVf^^, which occurs 
when the ST is at the edge of a harvesting zone, it follows that M = J^'-c , where {x] denotes the smallest 
integer larger than or equal to x. Note that M = 1 corresponds to the case where the battery is fully charged 
within one slot time; thus this case is referred to as single-slot charging. Similarly, the case of M = 2 is referred 
to as double-slot charging. It will be shown in this section that if M = 1 or M = 2, the battery power level can be 
exactly modeled by a finite-state Markov chain; hence, the transmission probability pt can be obtained. However, 
for multi-slot charging with M > 2, only upper and lower bounds on pt are obtained based on the Markov chain 
analysis for the case of M = 2. 

A. Single-Slot Charging (M — 1) 

If < P5 < rjPpVf^^, the battery of a ST is fully charged within a slot, i.e., M = 1. It thus follows that the 

battery power level can only be either or Pg at the beginning of each slot. Consider the finite-state Markov chain 

with state space {0, 1} with states and 1 denoting the battery level of power and P5, respectively. Furthermore, 

let Pi represent the state- transition probability matrix that can be obtained as 

^-Ph Ph 
Pg ^-Pg 

with Pg and ph given in (4) and (9), respectively. Then pt can be obtained by finding the steady-state probability 
of the assumed Markov chain, as given in the following proposition. 

Proposition 3.1: If < P5 < rjPprJ^^ or M = 1 (single-slot charging), the transmission probability of a typical 
ST is given by 



Pi = 



(12) 
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Proof: Let the steady- state probability of the two- state Markov chain be denoted by tti = [vri^o? ^1,1]? where 
TTi is the left eigenvector of Pi corresponding to the unit eigenvalue such that 

TTlPl =7ri. (15) 

From (15), the steady-state distribution of the battery power level at a typical ST is obtained as 

Pg Ph .... 

^1,0 = \ , TTl^l . (16) 

Vh ^Vg Ph+ Pg 

Note that the probability that a ST is fully charged at the beginning of each slot as defined in (10) is = tti^i 
in this case. Consequently, from (10), the desired result in (13) is obtained. ■ 
It is observed from (14) that in the single-slot charging case, pt depends only on Ap, and r^, but is not 
related to P5. The reason is that the battery of a ST is guaranteed to be fully charged over one slot if it gets into 
a harvesting zone; hence, the probability that a ST is fully charged = tti^i = does not depend on P5. 

B. Double-Slot Charging (M = 2) 

If TjPp^h^ < Ps ^ "^vPp^h^ or M = 2, a ST needs at most 2 slots of harvesting to make the battery fully 
charged. To establish the Markov chain model for this case, we divide the harvesting zone b(X^rh) into two 
disjoint regions, b{X^hi) and a(X, /ii, r/^), where hi = < '^h and a(T^x^y) = h{T^y)\h{T^x) denotes 

the annulus with radii Q < x < y centered at T G M^. It then follows that the region h{X^hi) consists of the 
locations at which the power harvested by a typical ST Y'^ from PT X is greater than or equal to Pg (i.e., single-slot 
charging region), while the region a(X, hi^r^) corresponds to the locations at which the power harvested by is 
greater than or equal to \Ps but smaller than Pg (see Fig. 3). For convenience, we define Hi = Uxe^p ^("^' ^1) 
and I-L2 = [Jxe^p ^(^5 hi^Vh). Note that H = T-Li UH2. We reasonably assume that Hi and H2 are disjoint since 
the harvesting zones are most likely disjoint as mentioned in Section II-B. 

Consider a 3-state Markov chain with state space {0, 1, 2}. Since the battery power level can only be either or 
in the range [^Ps, Ps] since rjPpVf^^ > in this case, we define state as the battery level of power 0, state 1 
with the power level in the range [^P^, P^), and state 2 with the power level equal to Pg. Note that in order to transit 
from state to 1, to 2, and 1 to 2, the harvested power at Y"" needs to be ^P^ < r]PpR~'^ < Ps, r]PpR~'^ > Ps, 
and rjPpR~^ > ^P^, respectively (or equivalently Y"^ needs to be inside H2, Hi, and H, respectively). Thanks to 
the fact that the minimum charging power is always larger than or equal to ^P^ in this case, we can determine 
the probability of the transition from state 1 to 2, i.e., from the battery power level in the range of [^Pg^Ps) to 
Ps, which occurs when Y^ is (anywhere) inside a harvesting zone (see Fig. 4(a)). Accordingly, the state- transition 
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Fig. 3. Divided harvesting zone for the case of double-slot charging (M = 2). 



probability matrix for the assumed 3-state Markov chain (see Fig. 4(b)) is given as 

1 - P/i P2 Pi 

l-Ph Ph , (17) 

Pg l-Pg _ 

where pi = Pr{y'^ G Hi} and p2 = Pr{y'^ G 7^2}. Notice that = P/i = 1 - e"^^'^^ since T^i U 7^2 = 

and we have assumed that %i and 7^2 are disjoint sets. Similarly to (7), the probability pi is given as 



= 1 - e~^^i^p, 



(18) 
(19) 



and p2 is given as 



P2=Ph- Pi 



(20) 
(21) 



Then we can obtain pt for this case as given in the following proposition. 

Proposition 3.2: If r]Ppr'^^ < Ps ^ '^vPp^h^ or M = 2 (double-slot charging), the transmission probability of 
a typical ST is given by 
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(a) Battery power state of ST (b) Markov chain model 

Fig. 4. The battery power state for the case of M = 2 and the corresponding 3-state Markov chain model, where (a) shows an example 
of the ST being in state 1 of the Markov model in (b), i.e., the current battery power level is in the range [^Ps,Ps)- 

Proof: The result in (22) can be obtained by following the similar procedure as in the proof of Proposition 3.1, 
i.e., by solving ^2^2 = 7r2, where 7r2 is the steady-state probability vector given by 7r2 = [7T2^o^ 7r2,i, 7r2,2]- Then, 
we obtain = 7r2,2 and then (22) is obtained from (10). ■ 
It is worth noting from (23) that pt in this case is a decreasing function of Pg since hi = (^^) " (23) is 
such a function. In other words, if Pg increases with fixed Pp and r/^, then the size of hi) (single-slot charging 
region) becomes smaller, which results in a ST harvesting for two slots to be fully charged more frequently, and 
thus a smaller pf. Hence, pt becomes smaller as well given pt = pfPg in (10). 

C. Multi-Slot Charging (M > 2) 

For multi-slot charging with Pg > 2irjPpr^^ or M > 2, the minimum charging power at the edge of the harvesting 
zone, ijPpT'^^, is smaller than ^Pg. Unlike the previous two cases of M = 1 and M = 2, the battery power level 
in this case cannot be characterized exactly by a finite-state Markov chain since it is not possible in general to 
uniquely determine the state-transition probabilities.^ However, we have shown that for the case of M = 2, the 
battery power level can indeed be characterized with a 3-state Markov chain regardless of the fact that we do not 
know the exact value of the battery power level in state 1, but rather only know its range [^Pg^Pg), provided that 
the minimum charging power rjPprJ^^ is no smaller than \Pg. Based on this result, we obtain both the upper and 
lower bounds on pt for the case with M > 2 as follows. 

As shown in Fig. 5, we divide the harvesting zone into 3 disjoint regions hi), a(X, /ii, /i2), and a(X, /i2, r/^), 
where < hi < h2 < Vh with hi given in the case of M = 2 and /i2 = (2^) " • ^^^^ ^^^^ ^("^' ^1) ^^^^ 

^For instance, if M = 3, following the previous two cases, we may divide the battery power level into 4 levels as 0, [|Ps, fi^s), 
[^Ps, Ps), and Ps and match each level to the states 0, 1, 2, and 3, respectively. Then it can be easily shown that the transition probabilities 
are unknown for some of the state transitions, e.g., from state 1 to 2. 
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Fig. 5. Divided harvesting zone for the case of M > 2. In this case, the amount of power harvested from PT X in a{X, h2,rh) is either 
overestimated as |Ps or underestimated as to obtain an upper/lower bound on pt in Section III-C. 



defined in tiie case of M = 2, wiiile tiie region a(X, /ii, /i2) consists of the locations in r/^) at which the power 
harvested from PT X is larger than or equal to ^Pg, but smaller than Pg, and the region a(X, /i2, rh) consists of the 
remaining locations in rh) at which the harvested power is smaller than ^Pg. Then, if we assume that the power 
harvested from a PT in the region a(X, /i2, ^/i) is equal to ^Pg (an overestimation), we can obtain an upper bound on 
Pt', however, if we assume it is equal to (an underestimation), we can then obtain a lower bound on pt, by applying 
a similar analysis over the 3-state Markov chain as for the case of M = 2. For convenience, we define the following 
mutually exclusive sets Ai = Uxe^p ^i)' ^2 = Uxe^p ^(^' ^i' ^2), and A3 = Uxe^p ^(^' ^2, ^/i), where 
Ai = Hi and AiU A2U A3 = H. Let p'^ = Pr{y^ e A2} and p3 = Pr{y^ G ^3}. It then follows that 
Pi P2 ~^ P3 ^ Ph^ where pi is given in (19) and 



= Pr{y^ G ^1 U A2} - Pr{y^ G ^1} = e'^^-^' - e'^^-^'^ 
P3=Ph-pi- P2 = e-"^-^^ - e-"^-^^ 

The following proposition is then obtained. 

Proposition 3.3: If Pg > 2r]Ppr'^^ or M > 2, the transmission probability of a ST is bounded as 

Pi + ^2 



(24) 
(25) 



<Pt < 



P2+P3 



Ph 



Pg- 



(26) 



{Pl +P2) +^41 + ^) P/z ( 1 + 

Proof: See Appendix A. ■ 

It is worth mentioning that the upper bound on pt is a decreasing function of Pg since hi = ^1^^ 

note that the bounds in (26) are tight in the case of M = 1 or M = 2, since P2 ^ Ps ^ with M = 1, and 
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Fig. 6. ST transmission probability pt versus ST transmission power Ps, with Xp = 0.01, rg — A, ru — 1-5, and Pp — 2. 

P2 = P2 and j>3 = with M = 2, thus leading to the same results in (13) and (22), respectively. 

Note that unlike the case of M = 2, it is not possible to verify in general whether pt for the case of M > 2 
is a decreasing function of Ps or not; however, it is conjectured to be so since a larger value of Ps will generally 
render a ST spend more time to be fully charged. We verify this by simulation in the following subsection (see 
Fig. 6). 

D. Numerical Example 

To verify the results on pt, we provide numerical examples as shown in Figs. 6, 7, and 8. For all of these 
examples, we set the path-loss exponent as a = 4 and the harvesting efficiency as 77 = 0.1. 

In Fig. 6, we show ST transmission probability pt versus ST transmission power Pg. It is observed that pt is 
constant if M = 1, but is a decreasing function of if M = 2, which agrees with the results in (14) and (23), 
respectively. It is also shown that if M > 2, pt is still a decreasing function of Pg as we conjectured. Moreover, 
the upper bound and lower bound on pt obtained in (26) for M > 2 are depicted in this figure. These bounds 
are observed to be tight when M = 1 and M = 2, while they get looser with increasing Pg when M > 2. The 
reason is that the size of the region a(X, /i2, Vh) shown in Fig. 5, in which we overestimate or underestimate the 
harvested power, enlarges with increasing Pg. However, since only small value of Pg is of our interest, we can 
assume that these bounds are reasonably accurate for small values of M. 
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Fig. 7. ST transmission probability pt versus PT density Xp, with rg — 3, rn — 1 and Pp — 1. 
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Fig. 8. ST transmission probability pt versus the radius of guard zone r^, with Xp — 0.01, rn — 1, and Pp — 1. 
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Fig. 7 shows pt versus PT density Xp, It is observed that for both M = 1 and M = 2, pt first increases with 
Xp when is small but starts to decrease with when Xp becomes sufficiently large. This can be explained as 
follows. If Xp is small, increasing Xp is more beneficial since each ST will get charged more frequently and thus 
be able to transmit (i.e., pf increases more substantially than the decrease of pg). However, after A^ exceeds a 
certain threshold, increasing Xp will more pronounce the effect of guard zones and thus make STs transmit less 
frequently (i.e., pg decreases more substantially than the increase of pf). 

In Fig. 8, we show pt versus the guard zone radius Vg. It is observed that pt is a decreasing function of Vg. 
Intuitively, this result is expected since larger Vg results in STs transmitting less frequently, i.e., smaller values of 
Pg, and it is known from (10) that pt = PfPg- 

IV. Outage Probability 

In this section, the outage probabilities of both the primary and secondary networks are studied. Let $t denote 
the point process of the active (transmitting) STs. In addition, let Ip and Is indicate the aggregate interference at 
the origin from all PTs and active STs, respectively, which are modeled by shot-noise processes [27], given by 
Ip = J2xe^p 9xPp\X\~^ and Is = J^Ye^t 9yPs\^\~^^ respectively. Note that in general, due to the presence of 
the guard zone and/or harvesting zone, in each time slot, the point process is not necessarily a HPPP; thus, 

(v) (s) 

Is is not the shot-noise process of a HPPP. Accordingly, the outage probabilities Po^t ^out primary and 
secondary networks, both related to J^, are difficult to be characterized exactly. To overcome this difficulty, we 
make the following assumption on the process of active STs, which will be verified later by simulations (see 
Fig. 10). 

Assumption 1: The point process of active STs $t is a HPPP with density ptXs. 

Let A(A) denote the HPPP with density A > 0. Under Assumption 1, the distribution of $t is the same as that 
of A{ptXs). It thus follows that Is can be rewritten as 

Is= J2 9yPs\Y\-''. (27) 

YeA{ptXs) 

(v) 

Consider first the outage probability of the primary network. Pout' which can be characterized by considering 
a typical PR located at the origin. Slivnyak's theorem [27] states that an additional PT corresponding to the PR 
at the origin does not affect the distribution of Therefore, the outage probability of the PR at the origin is 
expressed as 

Ip + Is + a^ 



Po-^^=Pr ,1 :,, <^P , (28) 
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Fig. 9. A typical SR located at the origin, for which there is no PT inside the shaded region b(Yo,rg). 

where gp is the channel power between the PR at the origin and its corresponding PT, and cr^ is the AWGN power. 

(v) 

Then, P^^;^ obtained in the following lemma. 

Lemma 4.1: Under Assumption 1, the outage probability of a typical PR at the origin is given by 

Pi^t^ = l-exp(-T^), (29) 

where 

rp-[\p+Pt\syyJ jOp4^+^^. (30) 

= 7rfr(^)r(l - I), with r(x) = J^y'^-^e-ydy denoting the Gamma function. 

Proof: See Appendix B. ■ 

(s) 

Next, consider the outage probability of the secondary network. Pout' which can be characterized by a typical 
SR located at the origin. Note that there must be an active ST, denoted by Yo, corresponding to the SR at the origin. 

(s) 

Since a ST cannot transmit if it is inside any guard zone, to accurately approximate P^^^ under Assumption 1, 
we consider the outage probability conditioned on that Yo is outside all the guard zones and thus there is no 
PT inside the disk of radius Tg centered at Yo (see Fig. 9). Let the event in the above condition be denoted by 
£ = {$p n 6(Yo, Tg) = 0}. Then the outage probability of a typical SR at the origin can be obtained as 



-a 
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Fig. 10. Outage probability of primary and secondary network versus SINR threshold, with a = 4, r] = 0.1, dp = ds = 0.5, = 3, 
rh = 1, Ap = 0.01, Xs 0.1, Pp 1, and Ps 0.1. 



where gs is the channel power between the SR at the origin and the corresponding ST Yq. From the law of total 
probability we have 



Pr{£} 



p(s, ^ L-^^--^" J ^ L^-;--'- _ . (32) 



Note that £ = n 6(1^05 7^ 0}- Then we have the following lemma. 

Lemma 4.2: Under Assumption 1, the outage probability of the typical SR at the origin is given by 

pis) ^ l-exp(-r.)-(l-p,) ^ ^33^ 
Pg 

where 

Ts = (^A, (gy^ +p,A,^ e|d2^+ ^£^. (34) 
Proof: See Appendix C. ■ 

(t?) (s) 

It is worth mentioning that the approximated result of Po^t ^out Lemmas 4.1 and 4.2, respectively, 
are valid only when Pp Ps, as assumed in this paper for the following reasons. First, to derive P^ut under 
Assumption 1, STs are uniformly located and thus can be inside the guard zone corresponding to the typical PR at 
the origin, and as a result cause interference to the PR. However, if we assume Pp » Pg, the interference due to 

(s) 

STs inside this guard zone is negligible and thus can be ignored. Next, to derive Pout' shown in Appendix C, 
the term Pr < / < 9s\£ \ in (32) can be assumed to be 1 only when Pp » P5. In Fig. 10, we compare the 
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1.2 




ST transmission power 

Fig. 11. Outage probability of primary and secondary network versus ST transmission power P^, with a = 4, r] = 0.1, dp = ds = 0.5, 

Tg =4,rh = 1, Xs = 0.1, Op = Os = 5, and Pp = 2. 

outage probabilities obtained by simulations and those based on the approximations in (29) and (33). It is observed 
that our approximations are quite accurate and thus Assumption 1 is validated. 

(s) 

In addition, it can be inferred from (33) and (34) that P^J^ is in general a decreasing function of Pg, since is a 
decreasing function of Pg. This implies that large ST transmission power Pg is beneficial to reducing the secondary 
network outage probability, although larger Pg also increases the interference level from other active STs. This can 
be explained by the fact that if Pg is increased, the increase of received signal power by the SR at the origin can 

be shown to be more significant than the increase of interference power from all other active STs. On the other 

(v) 

hand, from (29) and (30), it is analytically difficult to show whether P^^l is a decreasing or increasing function of 

Pg. This is because in general there is a trade-off for setting Pg to minimize the primary outage probability, since 

(v) 

larger Pg increases the interference level from active STs (resulting in larger Pout) but at the same time reduces 

(v) 

the ST transmission probability pt (see Fig. 6) and thus the number of active STs (resulting in smaller Pout)- 

(v) (s) 

Fig. 11, we show the outage probabilities Po^t and Po\jt versus Pg, with two different values of A^, respectively. 

(s) (v) 

It is observed that Po\jt is a decreasing function of Pg, whereas Pq^^ is quite insensitive to the change of Pg. 

V. Network Throughput Maximization 

In this section, the spatial throughput of the secondary network defined in (11) is investigated under the primary 
and secondary outage constraints. To be more specific, with fixed Pp, A^, r^, and r/^, the throughput of the secondary 
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network Cg is maximized over Pg and under given ep and 65. The optimization problem can thus be formulated 
as follows. 

(PI) : max. p,A, log2(l + (35) 
s.t. P^^<ep (36) 
P^t < es, (37) 

(v) (s) 

where Po^t P^J^ are given by (29) and (33), respectively. With other parameters being fixed, the transmission 
probability pt is in general a function of Pg (cf. Section III). Thus, we denote pt as pt{Ps) in the sequel. 

Since log2(l + Og) in (35) is a constant and P^ut' ^out monotonically increasing functions of Tp and r^, 
respectively (see (29) and (33)), (PI) is equivalently expressed as 

max. pt{Ps)Xs (38) 
s.t. Tp < fip (39) 
< /i^, (40) 

where /x^ = — ln(l — e^) and = ~ ^s)Pg)' Note that /ip and /i^ are increasing functions of and 

65, respectively. In general, it is challenging to find a closed-form solution for (38) with > 0. However, if we 
assume that the network is primarily interference-limited, by setting cr^ = 0, a closed-form solution for (PI) can 
be obtained as given in the following theorem. 

Theorem 5.1: Assuming = 0, the maximum throughput of the secondary network is given by 

where the optimal ST transmit power is 



log2(l + ^s), (41) 



Op y dp J y fJjp 

and the optimal ST density is 

_ fisif^p-^e^dlXp) 

Ag — 2 • l^-^J 

Pt{P:)es-djfip^ 

Proof: See Appendix D. ■ 
Note that since Pt{Ps) has been obtained in close-form for the case of < P* < 2riPpr^^ (i.e., M = 1 or M = 2 

in Section III), the optimal ST density A* in (43) can be obtained exactly for this case, according to (14) and (23). 

Otherwise, only upper and lower bounds on A* can be obtained, based on (26). 
Some remarks are in order. 
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Fig. 12. Optimal ST transmission power P* versus PT density Ap, with a = 4, dp = = 0.5, rn = 1, = 3, Pp = 2, = 0.3, and 
Op = Os = 5. 
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Fig. 13. Maximum secondary spatial throughput C* versus PT density Xp, with a = 4, dp = ds = 0.5, rn — l,rg — 3, Pp — 2, eg — 0.3, 
and Op — Os — 5. 
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14. Optimal ST density A* versus PT density Ap, with a = 4, dp = ds = 0.5, rh = 1, = 3, Pp = 2, eg = 0.3, and Op = Os = 5. 

It is worth noting that /i^ = — ln((l — es)pg) in (40) is an increasing function of PT density Ap, since pg given 
in (4) is a decreasing function of A^. Hence, the optimal ST transmission power P* given in (42) decreases 
with increasing Xp. This result is shown in Fig. 12, with three different values of Cp. 

In Fig. 13, we show the maximum secondary spatial throughput C| given in (41) versus A^ with Cp = 0.1, 
0.2, or 0.3. Note that from the perspective of RF energy harvesting, larger A^ is beneficial to the secondary 
network throughput. However, it is observed that if = 0.1, C| decreases with Xp, whlie for = 0.2 or 
0.3, C* first increases with Xp when A^ is small but eventually starts to decrease when Xp exceeds a certain 
threshold. The reason of this phenomenon can be explained as follows. When Cp is small as compared with 
Eg (e.g., = 0.1 in Fig. 13), the constraint in (39) prevails over that in (40), i.e., satisfying (39) is sufficient 
to satisfy (40), but not vice versa. Therefore, in this case, if Xp is increased, the active STs' density ptXg or 
Cg will be decreased to reduce Tp in (39), i.e., reducing the network interference level. However, when Cp is 
relatively larger (e.g., Cp = 0.2 or 0.3 in Fig. 13), (40) prevails over (39). As a result, if Xp is increased, then 
so is fig in (40), and thus ptXg or C| will be increased. However, if Xp exceeds a certain threshold, ptXg will 
be decreased to reduce in (40); as a result, C* decreases with increasing Xp. 

It is revealed from (43) that for given Xp, the optimal active STs' density pt(P*)A* is fixed under a given 
pair of primary and secondary outage constraints. In other words. A* is inversely proportional to pt{Pg). This 
implies that as pt converges to zero with A^ ^ (see Fig. 7), A* diverges to infinity at the same time, as 
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shown in Fig. 14. Thus, although the sparse PT density will lead to larger secondary network throughput (see 
Fig. 13), a correspondingly large number of STs need to be deployed to achieve the maximum throughput, 
each with a very small transmission probability As a result, only a small fraction of the STs could be active 
at any time, resulting in large delay for secondary transmissions or inefficient secondary network design. 

VI. Application and Extension 

In this section, we extend our results on the CR network to the application scenario depicted in Fig. 2, where a 
set of distributed wireless power chargers (WPCs) are deployed to power wireless information transmitters (WITs) 
in a sensor network. It is assumed that wireless power transmission from WPCs to WITs is over a dedicated band 
which is different from that for the information transfer, and thus does not interfere with wireless information 
receivers (WIRs). For simplicity, we assume that the path-loss exponents for both the power transmission and 
information transmission are equal to a. Moreover, the network models for WPCs and WITs as well as the energy 
harvesting and transmission models of WITs are similarly assumed as in Section II for PTs and STs in the CR 
setup. For convenience, we thus use the same symbol notations for PTs and STs to represent for WPCs and WITs, 
respectively. 

A. Transmission Probability 

Unlike the CR case, WITs in a sensor network do not need to be prevented from transmissions by guard zones, 
since there are no PTs present. As a result, a WIT can transmit at any time provided that it is fully charged. By 
letting Vg = 0, we have pg = 1, and from (14), (23) and (26) we obtain the transmission probability of a typical 
WIT in the following corollary. 

Corollary 6.1: The transmission probability of a typical WIT is given by 



1) If < < ryPpr-^ or M = 1, 



2) If ryPpr-^ <Ps< 2r]Ppr-'' or M = 2, 



Pt = T^. (44) 



Ph 
Ph + l + 



Pt = , , , p, - (45) 



Ph 

3) If > 2r]Ppr-'' or M > 2, 



P1+P2 ^ ^ Ph 

^-^*-rTTTS' ^^^^ 



where p/^ = 1 — e '^'^^K is given in (9); = 1 — e and P2 = e ^^l^p — e ^^1^^ are given in (19) and (21), 

respectively; P2 — e~^^p^^ — e~^^p^^ and p^ = e~^^p^^ — e~^^p^h are given in (24) and (25), respectively. 
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Fig. 15. WIT transmission probability pt versus WPC density Ap, with a — 4, r] — 0.1, rn — 1, Vg — 3, and Pp — 1. 

It is worth noting that unlike the CR setup, pt in this case is in general an increasing function of Ap since there 
are no guard zones and thus larger Ap always help charge WITs more frequently, as shown in Fig. 15. 

B. Network Throughput Maximization 

Note that unlike the CR setup, here we only need to consider the outage probability of a typical WIR at the 
origin due to the interference of other active WITs. Similar to Assumption 1, we assume that active WITs form a 
HPPP with density pt\s\ thus, the outage probability of a typical WIR at the origin can be obtained by simplifying 
Lemma 4.1 as 

i'iu1 = Pr|^^<^.l (47) 



where in this case is given by 



= l-exp(-r,), (48) 



r,^eldl^PtK + '^-^. (49) 

For the sensor network throughput maximization. Problem (PI) can be modified such that only the outage 
constraint for the WIR is applied. Thus we have the following simplified problem. 

(P2): max. p,(P,)A, log2(l + (50) 
s.t. P^,<es. (51) 
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The solution of (P2) is given in the following corollary, based on Theorem 5.1. 
Corollary 6.2: Assuming = 0, the maximum network throughput is given by 

C: = ^^log2(l + ^.), (52) 
where fi'g = — ln(l — eg), and the optimal solution (P*, A*) G M+ x R_|_ is any pair satisfying 

Ptip:)K = -4^- (53) 



^2 



Proof: With a"^ = 0, from (48) and (49), Problem (P2) can be equivalently rewritten as 



max pt{Ps)K (54) 

s.t. pt{Ps)\s<^^. (55) 

where /x'^ = — ln(l — 65). To maximize pt(Ps)Xs, then it is easy to see from (55) that the optimal solution is 
Pt{Ps)K = by multiplying it with log2(l + 6s), we then obtain C* in (52). 

■ 

Note that unlike the result in Theorem 5.1, the maximum network throughput remains constant regardless of 
Xp. This is because there is no primary outage constraint in this case and thus the optimal density of active WITs 
Pt(P*)A* is determined solely by the outage constraint of WIRs. On the other hand, if Xp is increased, we can 
effectively reduce the required WIT density A* for achieving the same C* since pt in general increases with Xp. 

VII. Conclusion 

In this paper, we proposed a novel network architecture enabling secondary users to harvest energy as well 
as reuse the spectrum of primary users in the CR network. Based on stochastic-geometry models and certain 
assumptions, our study revealed useful insights to optimally design the RF energy powered CR network. We 
derived the transmission probability of a secondary transmitter by considering the effects of both the guard zones 
and harvesting zones, and thereby characterized the maximum secondary network throughput under the given 
outage constrains for primary and secondary users, and the corresponding optimal secondary transmit power and 
transmitter density in closed-form. Moreover, we showed that our result can also be applied to the wireless sensor 
network powered by a distributed WPC network, or other similar wireless powered communication networks. 

Appendix A 
Proof of Proposition 3.3 

For both the upper and lower bounds, similar to the case of M = 2, we apply a 3-state Markov chain with state 
space {0, 1,2} with states 0, 1 and 2 denoting the battery power level of 0, in the range [^Ps^Ps), and equal to 
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Ps, respectively. First, consider the upper bound on pt. Since the harvested power in the region a(X, /i2,r/i) is 
assumed to be equal to ^Pg, it is easy to see that the state transition-probability matrix for this case is given by 



p(w) ^ 



l-Ph Ph 
Pg 1-pg 



(56) 



Let TT^^) = [ttq , tt}^ , ] denote the steady-state probability vector in this case. Solving 7r*^^)p*^^) = 7r^^\ 

(u) 

and thus the upper bound on pt can be obtained by multiplying with pg, 



we obtain tt. 



(u) 



Ph 



according to (10). 

Next, consider the lower bound on pt. Since the harvested power in the region a(X, /i2, rh) is assumed to be 0, 
it is easy to obtain the state transition-probability matrix for this case as 

l-(Pl+P2) P2 Pi 

P(')= 1-{P1+P'2) P1+P2 • (57) 

Pg ^-Pg _ 

Similarly to the derivation of the upper bound on pt, the lower bound on pt can be found by finding the 



corresponding steady-state probability ^ = 
of Proposition 3.3 is thus completed. 



P1+P2 



(Pl+P2)+P9[ 1 + 



- ^2 ^ 



, and then multiplying it with pg. The proof 



Appendix B 
Proof of Lemma 4.1 

(v) (v) 

For convenience, we derive the non-outage probability 1 — Pq^\. as follows with P^^t given in (28). 



pip) _ 

^out — 



gpPpdp ^ 



Ip + h + a^ - 
Pr{^P>^(/, + /. + <7^)} 



exp 



^p^p 2 
exp I 



(t , r I 

—B- [Ip + Is + (^ ) 



Opdp 
exp [ — -^Ip 



-Is 



exp - 



9 

JpUjrp 



(58) 
(59) 
(60) 
(61) 



where in (61), the expectations are separated since Ip and Is are assumed to be independent as a result of Assump- 



tion 1. Note that E/ 



exp 



and E/ 



exp 



are Laplace transforms in terms of the random 



variables Ip and Is, respectively, both with input parameter According to the result in [25, 3.21], the Laplace 

y rp 

transform of the shot-noise process of a HPPP A(A) with density A > 0, denoted by / = X^tgA(A) 
with input parameter s is given by 



E/[exp(-5/)] = exp(-(P5)' A(^), 



(62) 
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Fig. 16. Illustration of the optimal solution for Problem (PI). 



where {^r}TGA(A) is a set of i.i.d. exponential random variables with mean 1, and Lp is given in Lemma 4.1. Using 



(62), we can easily obtain E/ 



exp 



and E/ 



exp 



P - 



and by substituting them to (61), the 



proof of Lemma 4.1 is thus completed. 



Appendix C 
Proof of Lemma 4.2 



The term Pr 
given by 



Psd- 



< 



^1 in (32) is obtained by following the similar procedure in the proof of Lemma 4.1, 

9sPsd~ 



Pr 



< 



1 - exp(T5), 



(63) 



Ip + Is + (J^ 

where Tg is given in (34). 

Next, under the assumption Pp Pg, it is reasonable to assume that the interference from even only one 
single PT inside b{Yo,rg) is sufficient to cause an outage to the typical SR at the origin. Consequently, we have 
Pr I fj^jfl^2 < 6>5 1^} ^ 1. Substituting this result, (63) and Pt{£} = e"""^'^^ = - Pg into (32) yields (33). 
The proof of Lemma 4.2 is thus completed. 
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Appendix D 
Proof of Theorem 5.1 

From (30) and (34), the constraints < jip and < jig given in (39) and (40) are equivalent to pt{Ps)^s ^ 
h{Ps) and pt{Ps)K < f2{Ps). respectively, where 



h{Ps) = 



1 ^ , 



p 



As illustrated in Fig. 16, fi{Ps) decreases whereas f2{Ps) increases with growing P^. The shaded region in Fig. 16 
shows the admissible set of {Ps^ptXg) that satisfies the given outage probability constraints. It is observed that the 
optimal value of pt{Ps)Xs is the intersection of the two curves pt{Ps)Xs = fi{Ps) and pt{Ps)Xs = f2{Ps)' The 
intersection point can be found by solving fi{Ps) = f2{Ps)^ which has no closed-form solution in general with 
cr^ > 0. However, by letting = 0, the closed-form solution of P* can be obtained as ^ ^ ^^^^ 

2_ 2_ 

Pt{P^)Xl = /i(P;) and (64), we then obtain pt{Ps)K = and accordingly A* = ^^^^""^f 

Pt{P*)e^ dl^ipif 

Theorem 5.1 is thus proved. 
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